Bulletin of Electrical Engineering and Informatics 

Vol. 7, No. 4, December 2018, pp. 529-537 
ISSN: 2302-9285,001: 10.11591/eei.v7i4.1185 


□ 529 


Newly Developed Nonlinear Vehicle Model for an Active Anti¬ 
roll Bar System 


Noraishikin Zulkarnain^, Hairi Zamzuri^, Sarah ’Atifah Saruchi^, Mohd Marzuki Mustafa^, Siti 
Salasiah Mokri®, Nurbaiti Wahid®, Siti Nor Zawani Ahmmad’, Alias Jedi* 

'’“'’^Centre for Integrated Systems Engineering and Advanced Technologies (INTEGRA), 

Eaculty of Engineering and Built Environment, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia 
*’*Centre of Research in Engineering Education and Built Environment (PEKA) 

^Centre of Integrated Design for Advanced Mechanical Design (PRISMA, Eaculty of Engineering and Built 
Environment, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia 
^’^Vehicle System Engineering, Malaysia-Japan International Institute of Technology, Universiti Teknologi 

Malaysia, 54100 Kuala Lumpur, Malaysia 

’’Eaculty of Electrical Engineering, Universiti Teknologi MARA, Universiti Teknologi MARA, 23000 Dungun, 

Terengganu, Malaysia 

^Instrumentation & Control Engineering, Universiti Kuala Lumpur, Malaysian Institute of Industrial Technology, 81750 

Masai, Johor, Malaysia 


Article Info 


ABSTRACT 


Article history: 

Received March 30, 2018 
Revised Jul 24, 2018 
Accepted Nov 11, 2018 


Keywords: 

Anti-roll bar 
NSGA-II 
Optimization 
Ride and handling 
Vehicle model 


This paper presents the development of a newly developed nonlinear vehicle 
model is used in the validation process of the vehicle model. The parameters 
chosen in a newly developed vehicle model is developed based on CARSIM 
vehicle model by using non-dominated sorting genetic algorithm version II 
(NSGA-II) optimization method. The ride comfort and handling 
performances have been one of the main objective to fulfil the expectation of 
customers in the vehicle development. Full nonlinear vehicle model which 
consists of ride, handling and Magic tyre subsystems has been derived and 
developed in MATLAB/Simulink. Then, optimum values of the full 
nonlinear vehicle parameters are investigated by using NSGA-II. The two 
objective functions are established based on RMS error between simulation 
and benchmark system. A stiffen suspension provides good stability and 
handling during manoeuvres while softer suspension gives better ride quality. 
The final results indicated that the newly developed nonlinear vehicle model 
is behaving accurately with input ride and manoeuvre. The outputs trend are 
successfully replicated. 
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1. INTRODUCTION 

Vehicle suspension system is used to keep apart the vehicle from uncomfortable vibrations 
transmitted from the road excitation on the tyres and the driver. During a manoeuvre, it can keep the vehicle 
under control as it transmits the control forces back to the tyres [l]-[4]. The ride comfort and handling 
performances have been the major development objectives of vehicles in order to fulfil the expectation of 
customers in the development of mechanical and electronics vehicle technology. Thus, designing a suitable 
suspension system is always an important research issue for attaining the desired vehicle quality. Currently, 
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there are so many researchers in active suspension system field around the world that tackle these issues. 
However, some of them prefer to use active suspension others prefer to focus on the anti-roll bar itself 

Currently, an automotive engineer encounters many challenges in trying to have a good ride and 
handling whilst at the same time in cornering or straight line. Thus, this research proposed to investigate a 
new control design to provide the required active roll control of the anti-roll bar system for a ground vehicle. 
At the same time, to be able to find an improvement in handling capabilities without sacrificing the ride 
comfort. 

Due to the growing market share of light duty vehicle, automotive engineers have been prompted to 
examine body roll minimization strategies. During driving manoeuvres in high centre of gravity vehicles due 
to its roll-over, dangerous operating scenarios may induce drivers to drive in an aggressive way such as high 
lateral acceleration, rapid tire dilation and emergence lane change in [5]. Commercially, to counteract the roll 
movement, the most used topology is by using anti-roll bar. This bar can be implemented in a passive as well 
in active topology, where some of the most important parameters are again performance, energy consumption 
and costs [6]. 

From the previous researches, most researchers used genetic algorithm to optimise the system 
parameter such as routing in networks, job shop scheduling problem and isothermal liquid phase kinetic 
sequence [7]-[9]. Moreover, the Non-Dominated Sorting Genetic Algorithm II (NSGA-II) is proving to be a 
robust optimization algorithm for a wide range of multi-objective problems. The NSGA II Pareto ranking 
algorithm is an elitist system and maintains an external archive of the Pareto solutions [10]. 

Trade-off analysis is the one of the principles of active anti-roll bar system. It can help decision 
makers achieve the ride and handling performance requirements at the same time. However, the trade-off 
analysis in an active ant-roll bar system for tilt control design using NSGA-II are yet available in the 
literature [11]. Therefore, in this paper, a general methodology is proposed for conducting trade-off analysis 
for the selection problem of multi-objective of ride and handling for ground vehicle. A general trade-off 
based on multi-objective optimization framework for an active anti-roll bar system by optimising the novel 
control strategy parameters are established in this research [12]-[13]. 

Therefore, in this paper the newly developed nonlinear vehicle model for an active anti-roll bar 
system has been used NSGA-II optimization method for the future used of the controller design. This paper is 
organized as follows. In the next section, the full nonlinear vehicle model and vehicle parameters are 
presented. Then, the following section describes the newly developed nonlinear vehicle model. The fourth 
section shows the verification of newly developed vehicle model with CARSIM Model. Finally, the 
conclusion and future works are being concluded and suggested in the last section. 


2. FULL NONLINEAR VEHICLE MODEL 

A full nonlinear vehicle model described by using mathematical equations was developed as a 
MATLAB/Simulink model and was validated by comparing results from the double lane change test with 
speed 70 km/h. The steering input was used as the input to the simulation steering wheel angle and divided 
by a gain to represent the steer angle to the wheel in the vehicle model. The ride test was performed with 
CARSIM’s input which are roll bump signal, pitch bump signal and combination pitch and roll bump signal 
with 20 km/h. The whole vehicle model consists of five main sub-models i.e. ride model, tyre model, 
handling model, side slip angle model and longitudinal slip ratio model as illustrated in Figure 1. 



Figure 1. Full vehicle model 


3. NEWLY DEVELOPED NONLINEAR VEHICLE MODEL 

A new methodology is proposed to optimise the parameters of full nonlinear vehicle model based on 
ride and handling performances. The objective is to find the optimum values of the full nonlinear vehicle 
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parameters. The block diagram configuration of methodology on development of nonlinear vehicle model is 
illustrated in Eigure 2. 



Eigure 2. Newly developed vehicle model block diagram configuration 


The two objective functions are established based on RMS error between simulation and benchmark 
system. This is done by optimizing the nonlinear vehicle parameters by comparing with output response 
based on commercial available vehicle dynamic software (CARSIM). The vehicle dynamic software is used 
as a benchmark in order to provide similar results in vehicle model. Newly developed nonlinear vehicle 
model must fulfil some conflicting area. Among those is ride comfort which is attained by minimizing the 
error of body roll angle in speed bump test. While, a good handling of a vehicle is a desirable property which 
is achieved by minimizing the error of body roll angle in steering test. 

Objective 1 in the following form, which is consisted of RMS error value of the vehicle body roll 
angle for the handling performance means when the test input is handled in steering part as described in 
Equation 1. A small value of f-^ indicate good ride performance. Objective function 2 is built also using the 
mathematical formula of RMS error. RMS error value of the vehicle body roll angle for the ride performance 
means when the test input is handled in speed bump test as declared in Equation 3. A small value 
of /2 indicate good ride performance. 

The objective functions considered are minimization of (i) R.M.S error value of the vehicle body 
roll angle, which is for the handling performance on steering input test and (ii) R.M.S error value of the 
vehicle body roll angle, which is for the ride quality on the Speed Bump test. The following relations express 
these two objectives: 

(i.) Objl, fi=R.M.S error value of the vehicle body roll angle, which is for the handling performance on 
steering input test (Double Lane Change) 


f\ ~ l^rmsll 

where T is the transition time and the constraint function, 

f/i, A > 0 

Wi = •{ 

(0, Otherwise 


( 1 ) 

( 2 ) 


(ii.) Obj2,f2=R.M.S error value of the vehicle body roll angle, which is for the ride quality on the Speed 
Bump test 


ermsl = J^IoOf(t)dt - J^f^^0l(t)dt 
f 2 ~ I ®rms2 I 

where the constraint function is given by: 


(3) 

(4) 
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Tz, fi > 0 
. 0, Otherwise 


Therefore, the objective function with constraint is given by: 
F(x) ^ f^ + 


(5) 


Where is the penalty parameter which is for R-^ equal to 0.5 and R 2 also is 0.5 respectively. The value of 
penalty function is equal for both objectives in order to get equal trade-off between both objectives. In order 
to get the optimal solutions of parameter design for a newly developed nonlinear vehicle model, the RMS 
error of roll angle value is chosen as an objective function. From literature, the main goal of using anti-roll 
bar is to reduce the body roll [14]-[15]. For that reason, only the roll angle output is considered in 
optimization process. 

The newly developed nonlinear vehicle model be able to fine-tuned by choosing the suitable value 
of parameters are listed in Table 1. The eleven parameters chosen to optimise are roll axis moment inertia, /^, 
spring stiffness at the front right, K^fr, spring stiffness at the front left, K^fi, spring stiffness at the rear right, 
Ksrr, spring stiffness at the rear left, K^ri, damping stiffness at the front right, C^fr, damping stiffness at the 
front left, Cs/i^ damping stiffness at the rear right, damping stiffness at the rear left, roll spring 
stiffness, 1^0 and roll damping stiffness, C^. Table 1 shows the design variables, their bounds for newly 
developed nonlinear vehicle model and also initial design referring from work done by Hudha et. al. [16]. In 
this study, the value of 4 , K^fi, C^^i, C^. and K^. are observed as eleven 

design variables to be optimally found based on multi-objective optimization of two different objective 
functions [17]. The optimization problem is solved using Elitist Non-Dominated Sorting Genetic Algorithm 
(NSGA-II) [18]-[19]. Since it is a multi-objective optimization, a number of optimal solutions have been 
obtained and it is shown in the Figure 3 in which every solution is non-dominated by other solutions. 


Table 1. Design Variables and Lower-upper Bounds for Newly Developed Model 


Design Variable 

lower-upper 

Initial Design 

Ir 

100-1000 

700 

^sfr-~ ^sfh ^srr-~ ^srl 

500-1000 

750 

r r r r 

'-‘sfv’ ^s/1’ '-‘srv’ '-‘srl 

20000-40000 

30000 

C(f) 

100-4000 

3495.7 


100-60000 

56957 


These objective functions are considered in a Pareto optimization process to obtain some important 
trade-offs among the conflicting objectives, simultaneously. The evolutionary process of multi-objective 
optimization is accomplished with a population size of 30 which has been chosen with probability of 
crossover and probability of mutation as 0.8 and 0.1, respectively. These parameters setting of NSGA-II 
optimization method are shown in Table 2. A total number of 30 non-dominated optimum design points have 
been obtained. In order to analyse Pareto set and select good solutions, it is widely accepted that visualization 
tools are valuable to provide the decision maker in meaningful way. It is normally easy to make an accurate 
graphical analysis of the Pareto set points for a two-dimensional problem but for higher dimensions it 
becomes more difficult [20]-[22]. According to Pareto chart in Figure 3, a solution in central region of the 
trade-off is chosen for analysis. Therefore, the choice of parameters optimization for nonlinear model are 
based on the trade-off between both situations and gives the result in the Table 3. 


Table 2. NSGA-II User Defined Parameters 


Parameter setting 

Value 

Number of generation 

500 

Population Size 

30 

Probability of Crossover 

0.8 

Probability of Mutation 

0.1 

Distribution index in SBX 

20 

Distribution index in polynomial mutation 

20 
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Table 3. The Best Parameter Optimi zation of Newly Devel oped Nonlinear Model by NSGA-II 


Parameter 

Value 

Ir 

825 

^sfl 

798 

^srl 

802 

dsfl 

34560 

r 

'-‘srl 

36008 


57240 

^sfr 

829 

^srr 

813 

dsfr 

35019 

c 

'-‘srr 

35009 

C/j, 

3905 


Handling vs Ride using NSGA-II 



Handling Case, x 10”^ 


Eigure 3. RMS error roll angle in handling vs RMS error roll angle in ride: optimal front 


3.1. Verification of Newly Developed Vehicle Model with CARSIM Model 

The newly developed nonlinear vehicle model using NSGA-II optimization method is verified with 
CARSIMEd version 4.51, a well-known vehicle dynamics software in order to determine its output response 
effectiveness. The dynamics behaviour of a vehicle obtained from CARSIMEd. 

Is used as a benchmark for the 14 DOE nonlinear vehicle model under the same input conditions. 
Therefore, a performance comparison with benchmark vehicle is required. As a benchmark for validation 
purpose, a vehicle dynamic analysis software, CARSIM is used. The comparison between the full vehicle 
model with CARSIM benchmark is done by comparing the vehicle output response trend. 

This study is focused on both ride and handling parts of vehicle dynamic. Thus, the selected 
appropriate vehicle handling output responses are selected including roll angle and roll rate for validation is 
preferable. While, the chosen suitable vehicle ride output components also used roll angle and roll rate taken 
in consideration respectively. In the validation process, a comparison of the performance output trend of 
newly nonlinear vehicle model with CARSIM benchmark vehicle requires proper handling and ride test 
procedure. As such, two types of analysis procedures are used for the validation purpose which are roll mode 
bump and double lane change manoeuvre. 

3.1.1. Ride Analysis Result 

As for the ride test analysis, a rolling mode test is used [23]. Eor this, only one side of the vehicle 
hits the bump i.e. front left side and rear left side, to create rolling effect in the rolling mode test. The bumps 
were arranged in a way that it created the rolling effect. The dimension of the bumps used in this simulation 
is height at 0.075 m, weight at 0.11 m and length of 2.44 m. The vehicle was driven at a constant speed of 20 
km/h. Eigure 4 shows the illustrations of the test bumps. In the roll mode test, only the left side of the vehicle 
will hop on and hop down from the bump in order to model this condition, the front left road input Zy^i and 
the rear left road input Zyyi be given that the road input for the right hand side will be considered as zero. 
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Figure 5 shows the simulation and the benchmark results for this mode indicating good correlation in terms 
of similar trends. 



Figure 4. Rolling mode test 



(a) 

Roll bump test; Rolf rate 



(b) 

Figure 5. Roll bump test (a) roll angle response; (b) roll rate response 
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The body roll angle of Eigure 5 shows that the peak magnitudes from the full nonlinear vehicle 
model is almost identical with the CARSIM model response. The peak magnitudes at t =5.396 s is about 
0.1676 m/s^ and t =5.851 s is about 0.1564 m/s^ when the front left tyre hit the bump and the peak 
magnitudes at about 7.585 s is about 0.1529 m/s^when the rear left tyre hit the bump. The roll rate graph 
that is generated from the simulation of full nonlinear vehicle model with CARSIM model results show a 
good correlation. The data from 0 to 5 second can be ignored as the vehicle is accelerating in order to achieve 
a constant speed 20 km/h before hitting the bump. 

3.1.2. Handling Analysis Result 

Next, in the handling analysis, a double lane change manoeuvre test is used [23]. As for the double 
lane change test, it is taken from CARSIM where the steering wheel angle input for CARSIM model is 
exported to newly developed nonlinear vehicle model, thus both models were operated with the similar input 
condition. This manoeuvre test is operated in speed condition at 70 km/h. The full vehicle model 
performance with speed of 70 km/h is shown in Eigure 6. The newly developed nonlinear vehicle model 
outputs trend are successfully replicated. However, there are some error in magnitude for roll angle and roll 
rate. This error in magnitude occurs because of nonlinearity properties of tyre model. The results indicated 
that the newly developed nonlinear vehicle model is behaving accurately with input manoeuvre. 


Double Lane Change : Roll Angle 



(a) 


Double Lane Change : Roll rate 



(b) 

Eigure 6. Double lane change (a) roll angle response; (b roll rate response 
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4. CONCLUSION 

A newly developed nonlinear vehicle model for an active anti-roll bar system using NSGA-II 
optimization method is successfully presented in this paper. This methodology is proposed to optimise the 
parameters of full nonlinear vehicle model based on ride and handling performances. The newly developed 
nonlinear vehicle model using NSGA-II optimization method has been verified with CARSIMEd version 
4.51 that is a well-known vehicle dynamics software in order to determine its output response effectiveness. 
In the validation process, a comparison of the performance output trend of newly nonlinear vehicle model 
with CARSIM benchmark vehicle requires proper handling and ride test procedure. As such, two types of 
analysis procedures are used for the validation purpose which are roll mode bump and double lane change 
manoeuvre. The roll rate graph that is generated from the simulation of full nonlinear vehicle model with 
CARSIM model results show a good correlation in ride analysis part. However, there are some error in 
magnitude for roll angle and roll rate in handling analysis performance. This error in magnitude occurs 
because of nonlinearity properties of tyre model. For the conclusion, the results indicated that the newly 
developed nonlinear vehicle model is behaving accurately with input manoeuvre. 
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